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We study the lepton flavor violating (LFV) processes \i — ► &y, \i — ► 3e, and fi — > e 
conversion in nuclei in the left-right symmetric model without supersymmetry and 
perform the first complete computation of the LFV branching ratios B(fj> — > f) to 
■ leading non-trivial order in the ratio of left- and right-handed symmetry breaking 

scales. To this order, B(fi — > ery) and B([i — > e) are governed by the same combi- 
ng . nation of LFV violating couplings, and their ratio is naturally of order unity. We 
'i—i [ also find B(n — > 3e)/B(/j, — > e) ~ 100 under slightly stronger assumptions. Existing 
limits on the branching ratios already substantially constrain mass splittings and/or 
mixings in the heavy neutrino sector. When combined with future collider studies 
and precision electroweak measurements, improved limits on LFV processes will test 
the viability of low-scale, non-supersymmetric LFV scenarios. 
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Leptons of different flavors do not mix in the Standard Model (SM) of electroweak interac- 
tions as a consequence of vanishing neutrino masses. The observation of neutrino oscillations, 
however, has provided clear evidence that nature does not conserve lepton flavor and that 
the SM must be part of a more fundamental theory that allows for lepton flavor violation 
(LFV). In the most widely held theories of neutrino mass, LFV is generated at high scales 
that are well beyond the reach of present and future collider experiments. Searching for 
LFV among charged leptons in low-energy measurements is an alternative and important 
way to probe additional aspects of LFV at such high scales. 

Attempts to observe, and theoretically predict, the manifestations of LFV involving var- 
ious modes of muon decay have a long tradition. The rather small upper limit (2 x 10~ 5 ) on 
the branching ratio for the fi — > e + 7 decay determined by Lokanathan and Steinberger 
almost fifty years ago led to a flurry of theoretical activity (see e.g and many subsequent 
papers) that resulted in the realization that the electron and muon neutrinos are different 
particles - a fact confirmed experimentally shortly afterward. Over the intervening years 
the increase in the intensity of muon beams and advances of experimental techniques led 
to impressive improvement in the sensitivity of various searches for LFV. Even though no 
positive effects have been seen so far, the upper limits of the corresponding branching ratios 
became smaller by a factor of ~ 10 6 . At present, the most stringent limit on the branching 



ratio for \x — > cy is [3] 



= r r . (/ i + - e , +7 _\ < 1.2 x lCT 11 90%C.L. , (1) 

obtained by the MEGA collaboration, while for the process of /i — > e conversion in gold 
nuclei, the SINDRUM collaboration has obtained the limit 

r(/i- + A(N, Z)^e~ + A(N, Z)) 
'"- ~~ r(n-+A(Z,N)^v li + A{Z-l,N + i)) 
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The present limits on other branching ratios are similarly impressive: 1.0 x 10 12 for 
£ M +-e+ e -e+ @, 4.3 x 10~ 12 for 5^0, and 4.6 x lO" 11 for B^ e [7\. Two ambitious new ex- 
periments aiming at substantial improvement in the sensitivity are being developed: MEG 
plans to reach sensitivity of ~ 5 x 10~ 14 for -B„+_, e + +7 J8j, while MECO aims to reach 
~ 5 x 10~ 17 for B^ e in aluminum 

Theoretically, the focus in recent years has been on frameworks that could both account 
for neutrino mass generation at high scales and lead to observable LFV in future experiments 
with charged leptons. The direct effects of light neutrinos on charged lepton LFV are "GIM 
suppressed" by a factors of (Am^/M^) 2 ^ 10~ 50 in the rate and are, thus, entirely negligible. 
In order to obtain LFV effects that could be seen by experiment, a mechanism must exist 
for overcoming this GIM suppression. Such a mechanism necessarily involves physics at 
mass scales heavier than the weak scale. The primary motivation for LFV studies involving 
charged leptons is to help determine both the relevant scale as well as the most viable models 
associated with it. 

Although a variety of such models have been considered, based on various supersym- 
metry (SUSY) scenarios [HI 11, EH EH 13 EI|, or left-right symmetry [3, [I?}, the most 



commonly-quoted are SUSY grand unified theories (GUTs), wherein quarks and leptons are 
assigned to the same representation of the unification gauge group at the GUT scale. Conse- 
quently, the large Yukawa coupling responsible for the top quark mass also appears in LFV 
couplings [l3l. Eif . The latter then give rise - via renormalization group evolution - to sizable 
lepton flavor non-diagonal soft SUSY-breaking terms at the TeV scale. Superpartner loops 
that contain insertion of these terms then produce unsuppressed LFV transitions involving 
charged leptons. For example, in a SUSY SU(5) scenario, one has 13| 



fl™ = 2.4 x 10- 12 (3) 

neglecting gaugino masses. For superpartner masses of order the weak scale, one would 
expect to see non-zero signals in the up-coming MECO and MEG experiments under this 
scenario. Moreover, one would also expect to observe an order a suppression of B^ e relative 
to B^ ei in this case since the conversion process entails the exchange of a virtual gauge 
boson between leptons and the nucleus rather than emission of a real photon. 

In this paper, we study an alternative paradigm for LFV, wherein LFV occurs at much 
lower scales and does not require the presence of supersymmetric interactions to overcome the 
GIM suppression factor. In this scenario, neutrino mass generation occurs at the multi-TeV 
scale via a spontaneously broken extended gauge group 1 . LFV for charged leptons arises 
from the interactions of the additional gauge bosons, heavy neutrinos, and Higgs bosons 
associated with the extended gauge symmetry. As an explicit realization of this scenario, 
we work within the left-right symmetric model (LRSM) [20, 21, 2^1, which gives a minimal, 



non-supersymmetric extension of the SM with non-sterile, right-handed Majorana neutrinos. 
As such, it contains triplet Higgs fields that have non-zero hypercharge and that provide the 
simplest mechanism for generating a Majorana mass term 23[. As pointed out in Ref. 24]|, 



For a discussion of neutrino masses in this scenario and related phenomenological issues, see, e.g., Refs. 

HQ. 
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models of this type may give rise to unsuppressed operators for the LFV decay fi — > 3e, and 
these operators in turn induce logarithmically enhanced amplitudes for /i — > e conversion 
at loop level 2 . In the present case such effects - which were missed in earlier LRSM studies 



|16l Il7| - result from the presence of the triplet Higgs fields. The large logarithms can 
compensate for the O(a) suppression of B^ e relative to B^ ei that generically follows for 
SUSY GUTs, and for Higgs masses of order 10 TeV or below, both branching ratios may be 
large enough to be seen in future measurements. Roughly speaking, we find 



^,»ir'x W ^lHj' (5) 



}A1 _ m-7 „ „.|„ |2 



lTeV\V, Ml 



lO-'x^Pl^j (6) 




where M Wr is the mass of the right-handed charged gauge boson, M s ++ is the mass of a 
doubly charged, SU(2)r triplet Higgs, and 

9^ = Y,{^ R ) eN {K R ) N ^M N /M WR f , (7) 

with Kr being a flavor mixing matrix for the right-handed neutrinos of masses Mjv (see 
below). In the limit of degenerate, right handed neutrinos, the LFV factor = 0. For 
heavy masses at the TeV scale, present experimental limits already constrain this factor 
to be tiny: \g\ iv \ ^ 10 -2 . In this case, the heavy neutrino spectrum must either be nearly 
degenerate or devoid of significant flavor mixing. 

Note that both branching ratios are proportional to the same LFV factor, |<?if v | 2 . Naively, 
one would expect the loop graphs giving rise to /i — > (with heavy neutrino-gauge boson 
intermediate states) and the logarithmically enhanced loops that dominate /i — > e to have 
different prefactors. As discussed in more detail below, however, the doubly charged, triplet 
Higgs <5^ + and its left-handed companion can have LFV violating Yukawa couplings hij of 
(9(1), but the sum over intermediate states in the logarithmically enhanced loop graphs 
converts the sum over products of these couplings into g\f v . 

Should gib, turn out to be nonzero, then one would expect the two branching ratios to be 
of similar size since the product of the In 2 and a in B^ e is 0(1). We expect that any theory 
with non-sterile heavy Majorana neutrinos will contain such log enhancements, due to the 
presence of a more complicated Higgs sector than one finds in the SM. However, in SUSY 
GUT scenarios where LFV occurs at high scales, these logarithmically enhanced loop effects 
decouple below the GUT scale and do not affect the relative magnitudes of the branching 
ratios. Only when the symmetry-breaking scale is relatively light does one expect the two 
branching ratios to be commensurate in magnitude. 

Somewhat weaker statements about the relationship between B^ e and B^ 3e can also 
be made within the context of this model. In particular, we find 



\h„ P h 
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B^ 3e « 300 x x B^ e , (8) 



2 Specific realizations of these ideas been discussed for a doubly-charged scalar singlet |24| and R parity- 
violating SUSY [13. 
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so that if all of the triplet Higgs couplings hij are of roughly the same size and no cancellations 
occur in the sum g\ lY = J2j h^jh* e1 the \i — > 3e branching ratio should be roughly two orders 
of magnitude larger than the conversion ratio. Given the present experimental limits on 
■8(U-»3e> one would then expect B^ e to be of order 10~ 14 or smaller. As we discuss below, 
if the conversion ratio is found to be non-zero with significantly larger magnitude, then one 
would also expect to see a sizable effect in the channel r — > 3£ (where £ denotes a charged 
lepton). 

Finally, we observe that, while the logarithmic enhancement of Bz_^ e is a generic feature 
of any model that yields effective \i — > 3e operators at tree level, precise relationships 
between the various LFV observables depend on details of the model. In this respect, 
our perspective differs somewhat from the view in Ref. 0. Indeed, the presence of a 



common factor of |gif v | 2 in B^ n and - but not B^ 3e - and its relation to the 

heavy neutrino spectrum follows from the pattern of symmetry breaking in this scenario and 
the corresponding hierarchy of scales that enters the couplings of the right-handed gauge 
sector to matter. In order to implement this hierarchy in a self-consistent way, we adopt a 
power counting in k/vr, where Vr and k are the scales, respectively, at which SU(2) fi and 
electroweak symmetry are broken. In contrast to previous studies [3, [l7| , we compute all 
LFV contributions through leading, non-trivial order in k/vr and show that they decouple 
in the v r — > oo limit as one would expect on general grounds jjjj . In addition, we point 
out the prospective implications of other precision measurements and future collider studies 
for LFV in this scenario and vice-versa. The identification of such implications necessarily 
requires the adoption of a specific model, as the corresponding symmetries of the model 
dictate relationships between the coefficients of effective operators that would appear in an 
effective field theory framework. Thus, it is useful to have in hand a comprehensive treatment 
within various model frameworks in order to use experiment to discriminate among them. 
In R parity-violating SUSY, for example, the LFV couplings that generate \i — > e, etc. also 
appear, in general, in the mass matrices for light neutrino flavors [26], whereas in the LRSM 
LFV for charged leptons and light neutrinos are effectively independent. 

Our discussion of the calculation is organized in the remainder of the paper as follows. 
In Section II we review the main features of the LRSM and define the relevant quantities. 
In Section III the effective vertices are calculated and the effective Lagrangians for the LFV 
processes are determined. Some of the detailed formulae are collected in the Appendices. 
Section IV gives an analysis of the results, along with a discussion of the rates as well as 
their ratios. We conclude in Section V. 



II. THE MODEL 

The gauge group of the theory is SU(2) L x SU(2) R x U(1)b-l with the gauge couplings 
Ql = Qr = g for the two SU(2)s and g' for the U(l). In this paper we follow the notation 



developed in Ref. |27( where the LRSM, its quantization, and its Feynman rules are discussed 
in detail. Below, we give a very brief introduction to the model, and explicitly define the 
quantities used in subsequent analysis. 

The matter fields of the model include leptons (Ll,r) and quarks (Ql,r), which are placed 
in the following multiplets of the gauge group: 



1/2:0: -1), L lR = [ ,f ) : (0 : 1/2 : -1) 

R 
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QlL= \ d\) : (1/2 : : 1/3) ' QiR = [ d\ ) '■ ( ° : 1/2 : 1/3) ' (9) 

Here, i = 1,2,3 stands for generation number, and (II, Ir, Y = B — L) labels representation 
of the gauge group for each multiplet. The representation determines interactions of the 
multiplet with gauge fields. Before spontaneous symmetry breaking (SSB) the latter include 
(a = 1,2,3), and B» for SU(2) L , SU(2) R , and U(1) b _l gauge group factors, 

respectively. 

The SSB is achieved via the Higgs mechanism. The Higgs sector of the theory is not 
unique. However, the main results of this paper are largely independent of the details of 
the Higgs sector provided the LRSM has triplet Higgses and therefore heavy right-handed 
neutrinos. In our study we choose jU H2| a Higgs sector that consists of the bi-doublet 
: (1/2,1/2,0) and two triplets A L : (1,0,2) and A R : (0,1,2): 

(10) 

where the vacuum expectation values (VEVs) are shown in the second line. The most 
general Higgs potential with this field content has been analyzed in Ref. fl8l |. If one requires 
the scale Vr in the multi-TeV range (but not significantly larger), the only choice which 
avoids excessive fine-tuning and leads to acceptable phenomenology is to set to zero certain 
couplings in the Higgs potential as well as Vr llal . Moreover, we assume no explicit or 
spontaneous CP violation in the Higgs sector j28|. In summary, two distinct mass scales 
appear in the model: the electroweak symmetry breaking scale k ~ k% ~ «2 ~ 250 GeV, and 
the scale vr at which SU(2)r and U(1)b-l are spontaneously broken. Phenomenological 
considerations require vr ^> k. 




A. Physical fields 

After SSB matter and gauge fields acquire non-vanishing masses, which generally allow 
for mixing of the fields with the same quantum numbers. In the following, we identify masses 
and mixing angles which are important for our calculation. We omit the discussion of the 
quark sector of the model, as it is irrelevant for our work. 



1. Leptons 

The 3x3 mass matrix for charged leptons is Mi = {vd^2 + Vd^i) / v2, where y D and yjj 
are, respectively, the Yukawa coupling matrices for the bi-doublet <fi and its charge conjugate. 
Mi is diagonalized by a biunitary transformation Vj} M{V R = (Mi) diag . Since (Mi) diag k, 
one has yD,VD "C 1. Here, V[ R are 3x3 unitary matrices. These matrices relate the 
charged lepton mass eigenstates Il,r to the corresponding flavor eigenstates from Eq. Q: 
1'l,r - v l,rIl,r- 

Within the LRSM it is convenient to describe neutrino fields by four- component spinors 

, V%R = i<T*%L ■ (12) 




6 



The 6x6 mass matrix for the neutrinos is of see-saw type. It has both Majorana and Dirac 
entries: 

Mv = ( m\ Mr ) ' M ° = "71 <kVdKi + ^ DK ^ ' Mr = ^ VmVr ' ^ 
where yu is a 3 x 3 Majorana-type Yukawa coupling matrix. M v is diagonalized by a 
6x6 unitary matrix V: V T M U V = (M u )diag- This matrix relates neutrino mass and flavor 
eigenstates: n' L = V*Nl, n' R = VNr. Three eigenvalues of M v are small (denoted later by 
m u ), of the order of M^/Mr, and correspond to the light neutrinos of the SM. For Mr in the 
multi-TeV range, ide ntify ing Mp with the charged lepton mass clearly violates the 95% CL. 
limit J2 m is < 0.7 eV |29j from WMAP. However, appropriate choices of yc and yo (leading 
to Mp with non-trivial family structure and absolute scale on the order of m e j k ~ 10~ 6 ) 
result into m Vi consistent with present phenomenology. The remaining three eigenstates 
are predominantly right-handed with mass M n ~ M R ~ VmVr ~ 1 TeV, since we assume 
yu ~ C(l)- The amount of heavy-light mixing of the neutrino sector is set by the ratio 
e ~ Md/Mr ~ yD^/yMVR <C 1. As discussed below, we will expand our results in e and 
k/vr, and retain leading non- vanishing order. 

The LFV couplings of leptons to gauge and Higgs bosons are conveniently parameterized 
in terms of two 6x3 matrices |27( 



K L = V?VL, K R = V?V R , V=[ V f„\ . (14) 



At the leading order in e, the upper 3x3 block of Kl and the lower 3x3 block of Kr, 
respectively, describe flavor mixing in the light and heavy neutrino sectors. They are analo- 
gous to the CKM matrix which appears in the quark sector of the SM, and satisfy unitarity 
conditions up to corrections of order e 2 . In particular, the upper 3x3 block of is the 
familiar mixing matrix for light neutrinos [30j. As observed in the introduction, contribu- 
tions involving light neutrinos (and K£) to any LFV process are GIM suppressed relative to 
those involving heavy neutrinos. Therefore, the leading contributions to the LFV processes 
we consider depend on the masses and flavor mixing of heavy neutrinos only. 



2. Gauge fields 



The charged gauge bosons acquire the following mass matrix 



M 2 - — 



-2k\K2 

-2 , o„,2 



4 ^ -2kik 2 k% + 2v r ) ' ^ 

which is diagonalized via the mixing angle £ = — tan -1 (2kiK2/v r ) /2 with the eigenvalues 

+ Ak\ K^j /4. Here k + = yj n\ + re 2 , and the mass eigenstates 



g 2 ( «+ + v 2 R t 



R 



are related to the gauge eigenstates by 

W L = cos £ W± + sin£ W 2 , W R = - sin£ Wi + cos£ W 2 
The mass matrix for the neutral gauge bosons is 



M 2 



/ 9 -K 2 
2 ^+ 



2~ K + 







44 £ (4 + 4^) -2gg 



> v 2 



R 



-Zgg'vR 



(16) 



(17) 



2g' 2 v R J 
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It has the following non-zero eigenvalues (the third one is vanishing) 

\ W + + 2v R (g 2 + g' 2 ) 



Ml 



t V[s 2 4 + 2v Ua 2 + g' 2 )Y - W + 2<?'V+^ 



The explicit form of the unitary matrix that diagonalizes Mq is given in Ref . [27J . Here, we 
only list the expression for the Z\ — Zi mixing angle: 



1 . _i 5' 2 /tiv /cos 2&W 

— tjiri — 

2 2c? w {Ml-Ml) 



(18) 



(19) 



where 9w is the weak mixing angle, and cw(sw) is cos 9w (sin 9w) (we use this abbrevia- 
tion throughout). In this paper we work in the regime where 0, £ <C 1 3 , since both are 
0((k/v r ) 2 ). Note that because ys> <C yu we have e < (f>, £. 



3. Higgs fields 

With the Higg s sector described above, there are six neutral and four charged physical 
Higgs bosons [27j. At leading order in e, however, the neutral Higgs bosons do not contribute 
to the LFV processes involving charged leptons in the external states, and we do not consider 
them in the following. Two of the remaining bosons, Hf 2 -, are singly charged, with masses 
M Hl2 . The last two bosons, 6^^, are doubly charged, with masses M s ++ . Masses of the 

Higgs bosons depend on a number of parameters in the Higgs potential, with the natural 
scale M H ~ M s ~ v R [3, Ej. 



B. Lepton Interactions 

The LFV interactions of leptons with gauge (W 2 ), singly and doubly charged bosons are 
given by the following lagrangian densities: 



C 



cc 



N 



-fP R {K R ) 



l-W^ + l 



N-W. 



2 fi 



r — y 



9_ 

V2 



C s ±± 



9 
2 



Hf N (h P L ) I + H{ I (V P R j N 
6++ F (h L>R P LjR ) I + 5 L ; R I (h\ R P RjL ) V 



(20) 
(21) 
(22) 



where P L)R = (1 =F7 5 )/2, N = N L + N R = N c , I = Il + Ir and where we have neglected (D(£) 
terms. With our choice of the Higgs sector, it follows |l8( that the 3x3 matrix couplings 
Hl and h R can be identical (manifest left-right symmetry) or can have components differing 



3 The experimental limits on the gauge boson mixing angles are |£| < 3 x 10 3 and \<p\ < 1.8 x 10 3 
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by a sign (quasi- manifest left-right symmetry). In the manifest left-right symmetry case one 
finds 

M diag 

h L = h R = K T R -i—K R = h h = K* L h L . (23) 

M w 2 

Note that it is K R , not K' R that appears in the definition of h. Because Kr may contain 
Majorana phases, h is not necessarily proportional to the unit matrix even if all heavy 
neutrinos are degenerate. At leading order in e one has 



hij= E (#*) n< (^*L-V^n> (24) 

n=heavy 



"J 



tfh) = (~h r h) = £ x n (K R ) (K R ) =g lfv , (25) 

n=neavy 

where the sum is over the heavy neutrinos only. Eq. (|23j) relates the lepton-gauge boson 
couplings to the lepton-Higgs triplet couplings. We emphasize that it is specific to left-right 
symmetric models, and plays a central role in phenomenological applications. Generalization 
to the quasi-manifest left-right symmetry case is trivial and we have explicitly checked that 
a possible relative sign between Hl and tiR has no observable consequences in LFV processes. 
Finally, note that for degenerate heavy neutrinos, i.e., x n = const, one has gu v = 0(e 2 ) due 
to the approximate unitarity of the lower 3x3 block of Kr. Thus g^ v depends only on the 
mass square differences of the heavy neutrinos. The same is not true for the individual hijS. 



III. CALCULATION 



Within the LRSM we performed a complete calculation of the LFV muon processes 
fi — > ej, /i — > e conversion in nuclei, and /i — > 3e to leading order in the expansion parameters 
k/vr and e ~ M^/Mr. Diagrammatic contributions fall into three classes, schematically 
shown in Fig.0 Given the lagrangian in the physical fields basis [27], we have first identified 
all diagrams contributing to leading order in k/vr and Mo /Mr. We have then calculated the 
LFV vertices \i — > ej*, H eZ^, and /i — > eZ% [Fig. [T] (a)]. Finally, we have combined the 
LFV effective vertices with the gg7*, qqZ{ 2 interactions, and relevant box-type [Fig. U (b)] 
and tree-level diagrams [Fig. ^ (c)] to obtain the effective lagrangian for /i — > e conversion 
and n —>■ 3e. Our calculation and main results are described in this section, and some 
technical details are given in the appendixes El Gil 



A. Identifying the leading contributions 



In the LRSM at low energy all effects of the right-handed sector are suppressed by powers 
of k/vr, as a consequence of the decoupling theorem j25j. In our analysis we keep only the 
leading contributions in the expansion parameters k/vr and e. Throughout, we work in 't 
Hooft-Feynman gauge. Our findings can be summarized as follows: 
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FIG. 1: Diagrams contributing to /j, — ► cy, \i — ► e conversion in nuclei, and /i — > 3e. The wavy 
lines represent neutral gauge bosons (7 or Z\^)- fJ> — > e7 is described by class (a), /i — > e conversion 
is described by class (a) (attaching a quark line to the neutral gauge boson) and class (b) (with 
two external quark legs). /1 — > 3e receives in principle contributions from classes (a) (attaching an 
electron line to the gauge boson), (b), and (c). 




b) c) 

• The n — > ej* vertex receives its leading contributions at order e°, (k/vr) . In accor- 
dance with electromagnetic gauge invariance, however, the momentum-independent 
piece of the amplitude vanishes, so that the resulting vertex function is actually of 
order q 2 /v' R , q being the momentum transfer. Our expression is fully gauge invari- 
ant and respects the decoupling theorem. Note that when this vertex is inserted into 
the /1 — > e conversion amplitude, the q 2 in the vertex cancels the 1/q 2 of the photon 
propagator, leaving a contact interaction that scales as ~ 1/ v 2 R . 

• The \i — ► eZ| vertex is again not suppressed by powers of e or (k/vr), and we only 
keep the momentum-independent component. This result does not contradict the 
decoupling theorem, as one of the external states belongs to the heavy sector of the 
theory. Since q 2 « M§ , the contribution from this vertex to the \i — > e conversion 
amplitude goes as l/v R . 

• The n — ► eZ{ vertex nominally receives its leading contributions at order e°, (k/vr)°. 
However, the momentum-independent part of this class of diagrams sums to zero, 
in accordance with the decoupling theorem. We find that the leading non-vanishing 
contribution is 0((k/vr) 2 ). Consequently, the contribution from this vertex to the 
conversion amplitude is also ~ 1/^- 

For the kinematics of the LFV decays considered here, the momentum dependent contri- 
butions to the fi — > eZ{ 2 vertices are highly suppressed and can be neglected. 
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B. Effective vertices 



The fx — > e LFV vertices can be expressed in terms of known couplings and form factors 

F L,R, A L,R as follows, 

( Zi) = eG< _J_ , p (!) x 

V^Tt) 2 %C ff '"^ 1 ' RJ r 1 i 

v^tt) 2 s w c Wv ^oi2^ \ ^ R R)t*> 

-z8(47r) 2 m^g 1 '(A L P L + A i? P it )} / i , (29) 

where q = p e — Zfy, cr^ = fb^'Tf]- The /ie7* effective vertex has both "anapole" {F^ R ) 
and dipole (Al..r) terms 4 . Only the dipole terms contribute to the on-shell decay /J — > e7, 
while both anapole and dipole contribute to /x to e conversion in nuclei. 

The n — > e effective vertices receive contributions from the one-particle-irreducible dia- 
grams depicted in Fig. |21 as well as from external-leg corrections. The vertex corrections can 
be grouped into three classes: (i) gauge contributions (including unphysical Higgs exchange), 
(ii) singly charged physical Higgs contributions, and (iii) doubly charged Higgs contributions. 
Power counting implies that only certain combinations of gauge bosons, neutrinos, and Higgs 
particles contribute to leading order in k/vr and e. The relevant intermediate states are 
indicated diagram by diagram in Table ID 



TABLE I: Intermediate states contributing at leading order in k/vr and y£> to fJ-ej*, fieZl, 
and \i e Z^ effective vertices in 't Hooft-Feynman gauge. For each topology in Fig|2] we list the 
intermediate states as they appear starting from the muon vertex and following the loop counter- 
clockwise. Neutrinos are denoted by (heavy) and N\ (light). G\ t2 denote the unphysical Higgs 
fields associated with the longitudinal polarization of the gauge bosons W\ t2 . 





a) 


b) 


c) 


d) 


e) 


f) 


7 


W 2 ,W 2 ,N h 




W 2 ,G 2 ,N h 


G 2 ,W 2 ,N h 


G 2 ,G 2 , Nh 

HuHu iVj 
g±± g±± 


Hj H) U L,R 


Zi 


W 2 ,W 2 ,N h 


N h ,N h ,W 2 


W 2 ,G 2 ,N h 
W 2 ,H 2 ,N h 
W 2 ,G u N h 


G 2 ,W 2 ,N h 
H 2 ,W 2 ,N h 
Gt,W 2 ,N h 


G 2 , G 2 ,Nh 
H 2 ,H 2 ,Nh ; G 2 ,H 2 ,Nh ; H 2 ,G 2 ,N\ 1 
Gi,G\,Nh ; Gi,G 2 ,Nh ; G 2 ,G\,N\ l 


N h ,N h ,G 2 


z 2 


W 2 ,W 2 ,N h 


N h ,N h ,W 2 


W 2 ,G 2 ,N h 


G 2 ,W 2 ,N h 


G 2 ,G 2 , Nh 


Ni,Ni,G 2 



4 The first term in Eq. H29[) involves a coupling of the flavor- violating lepton current to the electromagnetic 
current rather than to a field associated with the corresponding vector potential. Zeldovich referred to 



this interaction as an anapole coupling 
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1. /j,eZ* vertex 

In this case the leading diagrams involving triplet Higgs (singly and doubly charged) sum 
to zero, and the main effect stems from gauge contributions. When working in 't Hooft- 
Feynman gauge, one needs to include the effect of unphysical Higgs exchange, and their 
mixing with other physical and unphysical scalars of the theory (terms proportional to 771 
and r)2 below). In terms of the heavy neutrino masses (M n ), mixing matrix Kr, and the 
ratios x n = (M n /M W2 ) 2 , y n = (M n /M Wl ) 2 , z n = (M n /M H2 ) 2 , the resulting form factors have 
the following structure: 



7(1) 

R 



,(1) 



n=heavy 



o 



rip 



Vo <Si (x n ) + 2 7] X D X (x n , y n ) + i] 2 D 1 ( 



(30) 
(31) 



where (k 2 _ 



K 2 K?,) 



sm 1 



v / cos~2l^ 

( «1«2 Y ^ 1 
\K + v R J ~~ 2 



M 2 W2 ' 




(32) 
(33) 



FIG. 2: Basic topologies for the one-particle-irreducible contributions to the fJ-ej*, fieZ*, and 
lieZ^ effective vertices (detailed version of Fig. U (a)). Wavy lines represent gauge bosons, dashed 
lines represent scalars (physical or unphysical), full lines represent leptons (charged or neutral). 
The internal particles contributing at leading order to each topology are listed in Table [IJ 
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K_ 



V2 
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Mi 



Hi 



and the functions S\(x), Di(x,y) are defined in Appendix 1X1 



(34) 



2. fi e Z| vertex 

As in the previous case, the leading term arises from gauge-lepton interactions, and the 
leading physical Higgs effects cancel out. With the notation established above, we find: 





— C W E 

n=heavy 
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' era \ / nfj, 



Si{x n ) 



(35) 
(36) 



3. fi e 7* vertex 



Both anapole and dipole transition form factors receive non-vanishing leading contribu- 
tion from gauge diagrams and exchange of singly and doubly charged triplet Higgs particles. 
Neglecting charged fermion masses (see Appendix El) and using Eq. (|25jh the various ampli- 
tudes read: 



A, 
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1 



up. 
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fi 



KwT 



E 

n=heavy 



en 



#1 



ra/i 



1< 

3M 2 ++ 



3 M 2 

°R 

24 M 2 + 

^2 1 



(37) 
(38) 
(39) 
(40) 



The functions 62,3(0;) are given explicitly in Appendix El The most important feature of 
these results is the logarithmic enhancement (q 2 ~ — m 2 ) of the anapole transition form 
factors, arising from the doubly charged triplet Higgs diagrams. This implies that in the 
left-right symmetry framework, \i — > e conversion in nuclei is as strong probe of LFV as 
fi — > cy since its amplitude is logarithmically enhanced, and thus compensates for the extra 
factor of ~ a. This effect was pointed out for a larger class of models in Ref. |24| within an 
effective field theory approach. Its consequences within the LRSM will be discussed in the 
next section in detail. 

In Appendix |B] we report full expressions for the ^e r j* form factors (including charged 
fermion masses) in terms of h, h (i.e. without using Eq. (|25jl). 
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a) b) 



FIG. 3: Box diagrams contributing to 4 J (Fig. □(b)). 



C. Effective Lagrangian for \i — > e conversion 



The effective Lagrangian for \i — > e conversion receives contributions from (i) tree level 
exchange of heavy neutral Higgs states; (ii) box diagrams depicted in Fig|31 (hi) LFV effective 
vertices, with the gauge boson attached to a quark line (the relevant quark-gauge couplings 
are summarized in Table ITT|) . Inspection of the neutral Higgs couplings implies that the 
ratio of effective couplings glf ee and gf^ generated by tree level Higgs exchange and loop 



2 /(a/An) -C 1. Therefore, we safely 



corrections, respectively, scales as gtree/ 9ioo P ~ (Vd) 
neglect the Yukawa suppressed tree level diagrams. 

After casting the fie Z± 2 vertices and the /i e 7* anapole terms in the form of a current- 
current interaction, the effective lagrangian can be written as in Ref. IIT 



AG F e 



^ea^(A L P L + A R P R )v-F^ 



Gi 



~^/2^ [9Lv(q)PL + 9Rv(q)PR 



+ejf, gLA(q)Ph + 9RA{q)PR n ® li^i^q \ + h.c. 



(41) 



where F^ u has to be understood as the classical field produced by the nucleus. In terms of 
the box contribution (S±(x) is defined in Appendix IA"|) 



? (B) 
R 



8 J2 (K R ) (k r ) S 4 (x n ) 

* — ' V / en \ / n/i 

n=heavy 



(42) 



and the LFV form factors, the couplings giv,Rv(q) ar e 

a 



9Lv(q) 
9Rv(q) 



4tt l q ' 



8ns 2 



w 



|-2sin 2 ^^) + l^ 



),(D 



4c 4 



w 



MKpiB) (B) 



(43) 



(44) 



The expressions for gLA,RA(q) are obtained by replacing i>W with a^' in gLV,Rv(q)- We remark 
that all the contributions to gRv{q) in Eq. f|44|) enter at leading order k/v r , contrary to what 
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TABLE II: Vector and Axial- Vector couplings of u and d quarks to Z\, and 7. We list for 
completeness the effective Vector and Axial- Vector couplings induced by box diagrams of Fig. [21 



Zi 


z 2 


7 


BOX 


„.W _ 1 _ 8 2 
u — J- 3 b W 

= 1 


- 1 - 

f« — 1 3*^/ 

flu = 1 "I - 


«u - 3 


= 1 

al B) = -1 


v d — 1 + 3 S W 

4 1} = -1 


„W - -1+4 2 
v d _ L ^ 3 S W 

(2) - I _ 2s 2 


V d - 3 


U d - 4 

a (B) - i 
a <2 _ 4 



appears in earlier calculations 0, E3|- I n R- e ^ |3| om y -^r 2 '' an -d -^r 8 ^ were included, while 
the authors of Ref. considered only Fj/ R . Both of these previous studies omitted the 

dominant, logarithmically-enhanced contributions from FjJ R . Finally, we note that upon 
taking matrix elements of C^ e in nuclei, the following combinations of gLV,Rv{q) become 
relevant: 

9lv,rv = 2gLV,nv(u) + 9LV,nv(d) , (45) 
9lv,rv = 9lv,rv(u) + 2g LV>RV (d) . (46) 



D. Effective Lagrangian for fi — > 3e 



The process fi — > 3e can occur in the LRSM through (i) tree level exchange of doubly 
charged Higgses (via the interaction of Eq.(j22J); (h) one-loop effective fi — > e vertex, with 
an electron line attached to the gauge boson; (iii) box diagrams. Barring the unnatural 
possibility that M s ±± 3> Mw 2 , the loop amplitudes (ii) and (iii) are suppressed by the 

L,R 

standard a/ir factor, and therefore in our analysis we disregard them. 

Doubly charged Higgs particles mediate at tree level also the decays r — > l a lbl c , with 
la,b,c — A 4 ) e - 111 compact notation, the effective lagrangian for four-lepton processes is given 
by: 



9' 



( ' D 



l iR IjR 



urI 



mR 



+ {L^R) 



(47) 



IV. ANALYSIS 



Based on the results described in the previous section, we now discuss the phenomenology 
of lepton flavor violation in muon decays within the LRSM. There are three main objectives 
of our analysis. First, we shall identify relations between LFV rates that are largely inde- 
pendent of the model parameters, and therefore can be considered as signatures of left-right 
symmetry broken at the multi-TeV scale. The pattern emerging is remarkably clear, and 
could be confronted with experimental findings in the next decade: the branching fractions 
for n — > e conversion and \x — > ej are expected to be very similar, and two order of mag- 
nitude smaller than the one for \x — > 3e (with some caveats). Second, we shall study the 
constraints on heavy neutrino masses and mixings implied by present experimental limits 
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on LFV processes. And third, we shall discuss the impact of future experiments, including 
collider measurements. 

Before describing the details of our analysis let us shortly recall the existing limits on 
the model parameters of interest to us. Direct searches imply that M W2 > 786 GeV, while 
singly- and doubly-charged Higgs particles should be heavier than ~100 GeV Indirect 
bounds are stronger and require the Higgs masses to be on the TeV scale. In summary, the 
existing phenomenology is consistent with the heavy sector masses being generically at the 
TeV scale or above. In what follows, we shall explore the consequences of a heavy mass scale 
being in the range 1-10 TeV, which can be tested in the foreseeable future. 



A. Setting the stage 

The quantities of primary interest to us are the branching ratios: 

iy-> ery) z r(/x -> 3e; 

-D^e 7 — (0) , -D^e — v z > ^V— 3e — ( 0) , ) 

1 fj, 1 capt 1 ^ 



where = (G^m^)/(1927r 3 ), and for the capture rate r^ pt we take the experimental 
values. The expression for the conversion rate T^ onv involves the overlap integrals 



= -±= r drr 2 N^pM (g^g- + f;f;) , (49) 



D 



2V2 Jo 

4m„ r°° 



jf drr 2 E(r) (//, / ;; • /, g~) . (50) 



Here = Z, N^ n ' = A — Z; p( p>n ) are proton and neutron densities, E[r) is the electric 
field generated by protons, and g~ f~, are the upper and lower components of the initial 
bound muon and final continuum electron wavefunctions, obtained by solving the Dirac 
equation. The overlap integrals have dimension of (massQf^ 2 , and in our study we use the 
numerical results for them reported in Table I of Ref. 32]. In terms of the form factors 
calculated above and D, V^ n \ the relevant branching fractions read: 

= 384vrV(|A L | 2 + \A R \ 2 ) , (51) 
B^ e = ^ (|A£D+S«V«+0My^ , (52 ) 

J- capt 

B - l \h h*\ 2 ( M ^ + M ^ \ (53) 

While B^se has a relatively simple structure, in general B jl ^ n and B^ e depend on a large 
number of unknown model parameters. However, under the rather natural assumption of 
a "commensurate mass spectrum" for the heavy sector of the model {i.e., My/ 2 ~ M s ++ ~ 
M s ++ ~ Mh+)i the problem becomes more tractable. Specifically, if M W2 , M s ++, M s ++, 
M H + , and the heavy neutrino masses M n are all of the same order of magnitude (in practice 
we shall assume 0.2 < Mi/Mj < 5 for each pair of masses), the amplitudes for fj, — > e and 
/i — > e7 become approximately proportional to gu v , defined in Eq. (|2*3jl . This is based on 
the following observations: 
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(i) Doubly charged Higgs contributions to the couplings A LR) gLvi.1) an d gRvio) are 
linear in x n (hence proportional to gif v ), and are sizable (the anapole transition form 
factor receives a large logarithmic enhancement). 

(ii) Gauge contributions depend on x n through the functions Si(x). These terms always 
represent a small correction to the Higgs contribution because 



m; 



(a) : \Si(x)\ « x £ log ^±±; (b) : S^x) « £ log 



3 m 



M] 



&++ . 



3 m 



2 



within the region (0.2) 2 < x < 3, where the lower limit follows from our assumption of 
commensurate spectrum and the upper limit from the vacuum stability condition [3^, 



Condition (a) ensures that gauge terms are small in the case of non-degenerate heavy neu- 
trinos, while condition (b) suppresses them in the case of nearly degenerate neutrinos. In 
what follows we account for the small gauge-induced contributions to the various couplings 
by expanding the Si(x) around x = 1.5, and keeping only the linear term. We have checked 
that the residual dependence on the expansion point x is small, and does not affect our 
discussion and results in a significant way. 

The above considerations about the relevance of g^ v remain true even in the unnatu- 
ral limit ^fj++ <S M W2 , but become invalid in the opposite limit M s ++^ 3> M W2 , as for 
M s ++ ~ 10M\y 2 the Higgs mass suppression compensates the logarithmic enhancement. 
Such unnatural limit will not be considered here. 

In summary, in the natural scenario of commensurate mass spectrum in the heavy sector, 
B^_> e and B^ ei are driven by a single combination of heavy neutrino masses and mixing 
parameters, which we defined as g^. Moreover, B^ e and B^ n depend only on four 
independent parameters (g\f v , M W2 , M s ++, M s ++), and have the generic structure 

L R 



d i ,2 ( M W2 M W2 _ M 5t+ _ M s ++ \ 

We shall next explore the consequences of such simplified form. 



B. fj, — > e conversion versus [i — > e^f 

The first important consequence is that the ratio R A = B^J ' B^ ei does not depend on 
Pifv, and is a function of \og(M W2 /M Wl ),r L ,r R . Our explicit analysis shows that, for input 
parameters in the commensurate range 0.2 < < 5, R A varies at most by 30% for any 
fixed value of M^/M^. As illustration, in Fig. 0]we show the ratio R for aluminum R Al as 
a function of My/ 2 /My/ 1 for a range of values (the variation with r R is much smaller). The 
most striking feature of our result is the near independence on the heavy mass parameters 
(as long as they stay in the natural range), leading to a distinctive prediction of the LRSM 
for R Al . This ratio is of 0(1) in this model and it is naturally confined between 1 and 
2, as shown by the gray area in FigEJ The absolute scale on this plot can be understood 



5 In the case of Al, the relevant conditions are | S3 (a;) | <§; x/3, and ^(a;)] <C 1/3. 
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as a consequence of the logarithmic enhancement of the anapole form factor contributing 
to B^ e . Different values of R A (in particular values smaller than unity) can be hardly 
accommodated without unnatural tuning of mass parameters. Indeed, for mass parameters 
just above the present direct limits (M W2 = 0.8 TeV and Mg++ = 200 GeV), we find 

R Al = 0.8, which can be considered the minimal acceptable value within this model. This 
prediction is substantially different from R-parity conserving SUSY scenarios, and can be 
hopefully tested by future measurements of B Al _^ e (MECO) and (MEG). 

The qualitative features encountered in the analysis of R Al apply to other elements as 
well. In particular, the ratio R A is always of 0(1). We have studied a few more examples, 
in the same range of mass parameters used above, finding: 

R Ti : 2 -> 3.5 , R Au : 2 -> 4 , R Pb : 1.5 -> 3 . (55) 

C. jjL — ► e conversion versus \i — ► 3e 

Under slightly stronger assumptions, it is also possible to derive an order-of-magnitude 
relation between B^ 3e and B^ e . Assuming dominance of logarithmic terms induced by 




o.5 ; 

100 200 300 400 500 

M W2 /M Wl 

FIG. 4: R Al = Bfi+JBn function of Mw 2 /M\y 1 , for different values of 77,. We keep 

r^j = 1, because the variation of i?" 4 ' with this parameter is considerably smaller then the variation 
with tl- The shaded band can be considered a prediction of left-right symmetry for R Al , assuming 
commensurate heavy sector. 
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Under the assumption that h^Jf^ 
between the three contributions, one then expects 



h^h*^ and that no cancellations occur 



B„ 



k f — I log 



9tt 2 



capt 



mi 



B 



»3e j 



(57) 



where fcy 
into 



\g\h\ 2 /\h^eh* e e\ 2 * s a number of order 1. For M 2 + 



fi— »3e 



3 x 10 2 



B 



Al 



ITeV, this translates 



(58) 



So, within this model, one expects that fi — > 3e could be the first rare muon decay to be 
observed. Sizable deviations from the above pattern would provide information about the 
parameters h^h^. In particular, B\^J B^ 3e ^> 1CT 2 would imply dominance of the I = fi 
and/or I = r contribution in \Ylihfj,ih* e \, and may lead to observable signals in r — > l a lbl c 
decays. On the other hand, B^ e / B^ 3e <C 1CT 3 would signal non-trivial relative phases 
among the couplings, necessary to suppress | J2i h^ih* le \ compared to \h^ e h* ee \. 



D. Constraints on heavy neutrino masses and mixing 

LFV in muon decays is driven by and the couplings related to heavy neutrino 
masses and mixing angles through Eqs. (|24|) - (|25|) . We now explore the correlations be- 
tween pifv and heavy mass parameters implied by present experimental limits and future 
limits/observations of £ff-> e > B At _ >e7 . Subsequently, we discuss the constraints on h^ e hl e 
implied by limits on B^ 3e . 



logio(fllfv) 
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FIG. 5: Correlations in the g\{ v -(M\Y 2 /M\Y 1 ) plane imposed by present and future (MEG and 

, represents the region allowed by limits 
Tr = 1. Lowering n and/or tr poses tighter 



MECO) limits on B^ e and B^ ei . The shaded area represents the region allowed by limits 
reported on the plot. In this plot we use tl 
constraints on gif v , for fixed Myy 2 . 
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In order to illustrate the generic model expectations for // — > e conversion and ji 
we show below approximate expressions for the rates (obtained by setting = tr = 



B 



1.5 x 1(T 7 |<7i fv | 2 



= X A x 10~ 7 



fi—te 



Star 



ZlTeVV 



ITeV 



Mi 



M 2 +4 

" 1 log f 

mi 



), 

(59) 
(60) 



where is the nucleus dependent numerical factor (we find = 0.8, 1.3, 1.6, and 1.1 for 
Al, Ti, Au, and Pb, respectively). These branching ratios have to be compared with present 
experimental limits: 



B 



<1.2xl0- n 3], Bll e < 4.3xl0- 12 [6], B^ e < 8 x 10~ 13 0, B^ e < 4.6 x lO" 11 . 

(61) 

Thus, assuming commensurate spectrum and g\{ v ~ 1 (i.e. large mixing angles and non- 
degenerate heavy neutrinos), consistency with present limits implies that the scale of SU(2)n 
breaking has to be around 20 TeV. On the other hand, for M\y 2 i n the 1-10 TeV range, present 
experimental limits already impose non-trivial constraints on g\f v (left panel in Fig.|SJ). Values 
of gif v at the 10~ 2 — 10 -3 level imply either small mixing angles in the heavy neutrino sector or 
nearly degenerate heavy neutrinos, on the scale set by M^ 2 . The most stringent constraints 
at present come from /i — > e conversion in gold. Future experiments MEG and MECO 
will be able to probe even higher mass scales and put more stringent upper limits on g\f v 
(right panel in Fig. |5J). Once again, /i — > e conversion will probe the model parameter space 
more strongly. 

Focusing on /i — > e conversion (present limits and projected MECO sensitivity), in Fig. EI 
we report a more detailed study of the constraints. At fixed Mw 2 , lowering or raising tl,r 
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FIG. 6: Correlations in the g\f v -(M\Y 2 /M\Y 1 ) plane imposed by B^ e , before and after MECO's 
goal has been reached, for different values of Tl = tr. The shaded area represents the region 
allowed by the assumed limits on B^f,. A non-zero mixing angle £ would further reduce the 
allowed region. As an illustration, the allowed region for £ = 10~ 3 is plotted in light-gray. 
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within the natural range 0.2 ^ n,R ~ 5, can change the bound on g lfv by an order of 
magnitude. Lighter Higgs particles imply tighter upper limits on gif v . Finally, the impact of 
a non-zero mixing angle £ (detectable, for example, through right-handed current signals in 
(3 decays) is also considered in FigJBJ A non-vanishing £ would imply [35| the upper bound 
Mw 2 /Mw 1 < 1/V£> an d thus narrow down the allowed region in the gu v - M W2 /M Wl plane 
(light-gray region in FigEJ). 

Additional information on heavy neutrino parameters can be obtained in principle from 
H — > 3e. The rate depends on doubly charged Higgs masses and the combination \h^ e hl e \ 
(Eq. 1(53)1 ) . The present limit B^ 3e < 10~ 12 implies (assuming M 5 ++ = M 5 ++) 6 

| Vft: j<1.55xlO-^(^) . (62) 

Thus, assuming M$++ ~ 1 TeV, the couplings ^ are constrained to be at the ~ 10~ 2 level. 
Unlike the case of g\{ v , however, the smallness of h^ e does not imply small mixing angles or 
almost-degenerate heavy neutrinos, because the Majorana phases contained in Kr may lead 
to cancellations in the sum of Eq. ([24)1 . 



E. Testing the model: interplay with collider measurements 



As noted above, information from LFV processes and other aspects of low energy phe- 
nomenology (such as signals of right-handed currents) can severely constrain the model 
parameter space in the near future. Moreover, given that B^ e and B^ n depend only on 
<7if v , M W2 , M s ++, M s ++, collider searches of heavy particles and low energy searches of LFV 
decays jointly provide a powerful probe of left-right symmetry. In fact, in the best-case sce- 
nario, separate measurements of B Al _> e , B^ ei and the mass parameters Mw 2 , M s ++ , M s ++ 
would allow one to test the model (four parameters versus five observables) . Even in less 
optimistic scenarios, one can imagine using collider information to narrow down the model 
predictions for LFV processes, or use observation of LFV to determine allowed regions in 
the heavy mass parameter space. 

As a simple illustration of this point, we show in Fig.[7|contour plots of R Al = B A L> e / B^ n 
in the M s ++-M s ++ plane, for two values of Mw 2 - We focus on the case of heavy masses in 
the 1-2 TeV range, which will be accessible at the LHC and Tevatron II [3^. In this 
mass-region, the model expectations are almost independent of My/ 2 . Moreover, one sees 
that values of R Al < 0.8 can only occur for < 100 GeV, already excluded by direct 

searches. Depending on future experimental developments, possible uses of the plots in 
Figs. [7| include: 

• Given measurements of Higgs and heavy gauge boson masses, one can infer rather 
precisely where to expect R Al within this scenario. 



6 A weaker upper limit on the same combination of parameters can be derived from searches of muonium 
anti-muonium transition In general, present limits on the flavor diagonal coupling h ee from 

Bhabha scattering 38], and other combinations of hij from rare r decays are much weaker (typically 

b t ^i Mc < ltr 6 
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FIG. 7: Contour plot of R Al = B^_ e /B^ ei in the M s ++ -M s ++ plane, for M W2 = 1 TeV (left 

panel) and for M\y 2 = 5 TeV (right panel). Each curve is labeled by the corresponding R Al . As a 
function of the Higgs mass along the line M x ++ = M x ++ , R Al reaches a maximum at M x ++ ~ 2 Mjy, 

°L °R °L 

and then decreases, due to decoupling of doubly charged Higgs bosons (the latter effect is not visible 
in the plots). 



Given an experimental signal for B Al _^ e and B^ ei , one can identify the allowed region 
in the M s ++-M s ++ plane, for different values of Mw 2 - Collider searches could then 
confirm or falsify the model expectations. As can be seen from the plots, however, 
in order to have a significant test, the fractional uncertainty on R Al should be at 
most 20% (otherwise most of the M s ++-M s ++ would be allowed). Given the projected 
sensitivities, this may be achieved at the next generation experiments if > 
2. x lCT 13 . 



V. CONCLUSIONS 

The study of flavor violation among leptons now lies at the forefront of particle and nuclear 
physics. The tiny masses of the three lightest neutrinos and the nearly maximal mixing 
among them stands in stark contrast with the situation involving quarks, and the origin 
of this difference remains a fundamental and unsolved puzzle. A variety of scenarios have 
been proposed that attempt to answer this question, and these ideas would have predictable 
consequences for other observables. In this study, we have analyzed the consequences of 
one such scenario - the left-right symmetric model - that entails a minimal extension of 
the SM gauge symmetries and that includes non-sterile, right-handed neutrinos whose mass 
could be generated at the multi-TeV scale, albeit with some fine-tuning. We have shown 
how it implies relationships among various LFV decays of the muon that could distinguish it 
experimentally from other models of LFV. We have also illustrated how direct searches for 
right-handed gauge bosons and triplet Higgs at the Tevatron and LHC would complement 
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the charged lepton LFV studies and either help favor or rule out the possibility of rather 
low-scale LFV without SUSY. 

The main conclusions of our study are: 

• The branching ratios B^_ e and B^ ei are similar in magnitude, in distinction to other 
possible scenarios which predict that B^ e /B^ ei ~ a. 

• Within the LRSM, and with reasonable additional assumptions, B^ 3e / B^ e ~ 300, 
making the process /i — > 3e perhaps easiest to observe. 

• The existing limits on the LFV muon decays already substantially constrain the mixing 
and mass splittings of the heavy right-handed neutrinos. The planned more sensitive 
experiments will therefore test the LRSM severely. 

If the LRSM scenario turns out to be correct, the deeper connections between the heavy 
and light neutrino spectrum would, then, have to be pursued by additional experimental 
and theoretical work. On the other hand, should experiment eliminate the possibility of 
non-supersymmetric, low-scale LFV based on the considerations discussed above, the lep- 
ton flavor problem will nevertheless remain a rich area of study, both theoretically and 
experimentally, for some time to come. 
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APPENDIX A: LOOP FUNCTIONS 



We collect here the functions Si(x) and Di(x,y) appearing in the expression of various 
ji — > e form factors. 



Si(x) 
S 2 (x) 

Ss(x) 
S 4 (x) 



Ax 



6 - 7x + x 1 + (2 + 3x) log x 



(1-x) 2 

x(4-3x) 2x(12 - lOx + x 2 ) 
(1-x) 2 3(1 - a;) 2 

x(l + 2x) 3x 2 
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Both Si(x) and Di(x, y) are regular at x = 1 and y = 1. Note that the potentially dangerous 
contribution involving the large mass-ratio i/ n = (M n /M^) 2 has a finite limit for y n — > oo: 

lim D\(x,y) = —7S&{x) . 



APPENDIX B: FULL EXPRESSIONS FOR ^7* FORM FACTORS 



In terms of the interactions vertices reported in Eq. (|21| l2*2*|) , and without neglecting the 
charged lepton mass-dependence of loops, the photonic form factors read: 
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Mk S 2 (x n ) 
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16 S B 
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4 + 



-16 5. 



mt 
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16n 2 A L 
16n 2 A R 
where the function S 5 (x) is: 
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24 V y 



^ M 2 + 
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S 5 {x) = J o dyy(l-y) log (x + y(l - y) 



(Bl) 
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(B3) 
(B4) 



(B5) 
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